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Abstract 
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^^ \ The occurrence of chaos for test particles moving around a slowly rotating 

black hole with a dipolar halo is studied using Poincare sections. We find a 

^ \ novel effect, particles with angular momentum opposite to the black hole 

(<— ^ \ rotation have larger chaotic regions in phase space than particles initially 

(n' 

moving in the same direction. 
C\ . 04.20.Jb, 04.70.Bw, 05.45.+b. 

bJO. 
> 



Typeset using REVT^ 



*e-mail: letelier@ime.unicamp.br 
''^e-mail: vieira@ime.unicamp.br 



The main paradigm in the study of the motion of stars in a galaxy is a model of a central 
bulge surrounded by a halo [|T|. In this context, for the particular axially symmetric case, 
arose the celebrated Henon-Heiles model [0] whose study has been the source of inspiration 
of many researches on chaotic behavior . The underlying theory in this case is the usual 
Newtonian Gravitation that for large masses and velocities is known to be less appropriate 
than Einsteinian General Relativity. In the late case the Newtonian potential is replaced by 
the spacetime metric and Newton motion equations by geodesies. This change of dynamics 
can produce dramatic effects, for instance, test particles moving in the presence of systems of 
masses that are integrable in Newtonian theory are chaotic in General Relativity, examples 
are: the fixed two body problem @J^, and particles moving in a monopolar center of attrac- 
tion surrounded by a dipolar halo 0. Also gravitational waves, a non existing phenomenon 
in the Newtonian realm, can produce irregular motion of test particles orbiting around a 
static black hole |0J^. Another distinctive feature of general relativity is the dragging of 
inertial frames due to mass rotation. This fact is observed, for instance, in the impressive 
differences of the geodesic motion in Schwarzschild and Kerr geometries |^. 

In this Letter we study the effects of rotation in the motion of a particle orbitating 
around a slowly rotating black hole surrounded by a dipolar halo. The non rotating case 
is chaotic, so we shall study mainly the change in the chaotic behavior due to the rotation 
of the center of attraction. A typical situation is represented by a galaxy with a rapid 
rotating center surrounded by a distant massive halo, ring or other shell-like distributions 
of matter. We study the motion of particles moving between the center and the halo whose 
first contribution is dipolar. This contribution is always present whenever the halo does not 
possess reflection symmetry with respect to the black hole equatorial plane. 

The metric that represents the superposition of a Kerr black hole and a dipole along the 
rotation axis is a stationary axially symmetric spacetime. The vacuum Einstein equations 
for this class of spacetimes is an integrable system of equations that is closely related to 



the principal sigma model ||TD[. Techniques to actually find the solutions are Backlund 



transformations and the inverse scattering method, also a third method constructed with 



elements of the previous two is the "vesture method" , all these methods are closely related 
p!0| . The general metric that represents the nonlinear superposition of a Kerr solution with a 
Weyl solution, in particular, with a multipolar expansion can be found by using the "inverse 
scattering method" [O. We find 



ds^ = gttir, z)df + 2gt^{r, z)dtd(j) + g^4,{r, z)d4P' + /(r, z){dz^ + (ir^), (1) 

where 
gtt = -[e~^Dpuv^^-iDpu^^-%Dpv^^ ^ ^Yiu" - 1) + e-«^P(p - ifiu" - 1) 

gt^ = 2pge-2^P("+"+^)[e-^^P"(e-^^P''(p + l)(n - v){u + 1) 

+e-4^P(p -!)(« + v){u - l)){v + 1)) 
-e-^Dpv^^ + 1)(m + t;)(M + l)(t; - 1) - e-^^P(p -l){u- v){u -l){v- 1)]/F, 
^w = [i.9tvf - {T^rf]l9tt, 

f = ^2^^^Dpi(v+l)u-v+l)^^-4Dpuv^p _ ^)2(^ _ ^)2 _^ g-4Dp({.+l)-a+.-l) (^ ^ ^)2^2] 

_^^-4.Dp({v+l)u+v-l)~iDpuv ^ _^ -|^\2/ _^ -|^n2 
_2Q-'^Dp{{v+l)u+v-l)-2Dp{iv+l)u-v+l)-4.Dpuv^^2 _ ^2n 2 _^ ^-SBp^^/ _ if n^) I H 

F = e-^^^"[e-«^^^(p + l)\u + 1)2 + e-«^P(p - 1)2(m - 1)^ - 2e-'^^''P^^+\u' - v^)q% 
H = 4exp[-2i)p((f + l)u + v-l)- 2Dp{{v + l)u - v + 1) 

-6Dpuv-{u^-l){v^ -1)DY]- (2) 

The coordinates (r, (p, z) are dimensionless and have the range of the usual cylindrical co- 
ordinates. They are related to u and v by: z = uv and r = [v? — 1)^/^(1 — v'^Y^'^ ^ u > 1 
and — 1 < f < — 1. Our units are such that c = G = 1; D represents the dipole strength 
and a the rotation parameter, q = a/m and p^ + q^ = 1. The coordinate transformation 
t' = t + 2aip, u = R/m — l,v = cos^, ^' = ^ reduces (0) with D = to the Kerr solution in 
the usual Boyer-Lindquist coordinates |]T2|. 



To study the slow rotation case is better to use the metric obtained by keeping the first 
order terms in the rotation parameter a in the exact metric (0). This approximation, for the 
parameters and range of coordinates used, will not produce a significant information loss; 
we shall comeback to this point later. We find for g^^ = g'^^ + ag^^^, 

u — 1 r-, \ 

9tt = —-exp{-2Duv), 

u + I 

9t<t> = ^t[(" + ^)(1 - ^) exp(-2L'(l -u- v)), 
+ {u -v){l + v) exp(-2Z)(l + u- v))], 

g<j}<j) = m'^{l - f^)(l + uy exp{2Duv), 

ill + T|2 
/ = m=^ ^" ; exp(i)[(M2 - l)(t;2 - l)D + 2uv - Av + 4]). (3) 

u — V 

The geodesic equations for the metric (|l]) can be cast as 

i = g'^Eb, = g'^'Eb, (4) 

f = -^XgfEaE, + /,,(r2 _ .2) ^ 2/,,fi], (5) 

^ = -^[^?^«^^ + fA^^ - ^') + 2/,rri], (6) 

where the dots denote derivation with respect to s and the indices a and h take the values 
it,(f)), g"'^ stands for the inverse of gab- Et = —E and E^ = L are integration constants; 
E and L are the test particle energy and angular momentum, respectively. The set (^-(§|) 
admits a third integration constant 

E3 = g'''EaEb + fir' + z') = -l. (7) 

Thus to have complete integrability we need one more independent constant of integration. 
In the case of pure Kerr solution [D = 0) we have a fourth constant due to the existence of 
a Killing tensor and for the non rotating case we have another constant related to a third 
Killing vector associated to spherical symmetry 0. 

The system d^)-® can be written as a four dimensional dynamical system in the vari- 
ables (r, z, Pr = r,Pz = z). A convenient method to study qualitative aspects of this system 



is to compute the Poincare sections through the plane z = 0. The intersection of the orbits 
with this plane will be studied in some detail for bounded motions. We shall numerically 
solve the system (§)-(|^) and use the integral (|^ to control the accumulated error along the 
integration; we shall return to this point later. 

The Poincare section for different initial conditions with energy E = 0.965 and angular 
momentum L = —3.75 (counter rotation) moving in an approximate Kerr geometry {D = 
0) with rotation parameter a = 0.01 and mass m = 1 (this value for the mass will be 
kept unchanged from now on) are presented in Fig. 1. We have the typical section of an 
integrable motion, i.e., the sectioning of invariant tori, for integrability and KAM theory. 



see for instance flS 



We also studied the same case for direct rotation L = 3.75, as well as, the corresponding 
Schwarzschild limit a = 0, and L = —3.75. All these cases present Poincare sections almost 
identical to Fig. 1. The section area for the Schwarzschild case is slightly smaller than the 
section area of Fig. 1 and the section area for particles in direct rotation in a Kerr geometry 
is even smaller. We have that counter rotation enlarges the area of the section and direct 
rotation shrinks it. This is clearly an effect of the dragging of inertial frames due to rotation. 

The motion of test particles around a static black hole with a dipolar halo (a = and 
D ^ ) is chaotic and it was studied in some detail in ^ for a different energy shell. In 
Fig. 2 we show the Poincare section for D = 0.0005 and the same values oi E = 0.965 and 
L = 3.75 as in Fig. 1. We find islands of integrability surrounded by chaotic motion. The 
two isolated islands around the points (10, 0.05) and (5, 0.075) are parts of the same torus. 
In the case studied in they were closer. 

Now we shall consider a particle moving around an slowly rotating attractive center 
surrounded by a dipolar halo for both direct rotation and counter rotation. In Fig. 3 we 
draw the Poincare section for D = 0.0005, a = 0.01, E = 0.965 and L = 3.75 (direct 
rotation). We see that the islands of stability are larger in this case than in the non rotating 
case (see Fig. 2); also we have new systems of small islands immersed in the chaotic region. 
We have that the chaotic region is smaller in this case than in the equivalent non rotating 



one. It does look like that the direct rotation diminishes the effect of the dipolar strength 
as a chaos source. In Fig. 4 we present the section with the same parameters of Fig. 3, 
except that now we have counter rotation, L = —3.75. In Fig. 4 we observe that the chaotic 
region increases in a significant way and also that the islands located around the points (10, 
0.05) and (5, 0.075) in Fig. 2 and 3 has disappeared in this scale. In other words, with 
D ^ 0, the counter rotating motion of particles is more chaotic than the static case, the 
later being more chaotic than the direct rotation case. Also, we can say that the counter 
(direct) rotation reinforces (weakens) the strength of the dipole as a source of the chaotic 
motion (for D = we have an integrable system). This effect is a manifestation of the fact 
that particles moving in non equatorial orbits of a rotating central body can suffer repulsive 
forces due to rotation. This fact allows in a Kerr geometry closed orbits of test particles 



moving on planes parallel to the equatorial plane [Q, though these orbits are not stable. 
Although in the present Letter we present results for particular values of the parameters 
involved, we did a rather extended numerical study that supports our conclusions. Figs. 3, 
and 4 being representative of this search. 

For the values of the parameters D = 0.0005, a = 0.01, E = 0.965 and L = ±3.75, we 
have that the particles move in the "box" 4.7 < r < 21, —5 < z < 9. We take as a measure 
of error the quantities 

^9a, = \{9:i-9a,)/g:u, A/ = \{r - f)/r\, (») 

where in these expressions the sum rule of repeated indices does not apply, g'^l and g^^, 
refer to the solutions (H) and (||), respectively. We find that for the above mentioned range 
of coordinates and the values of parameters used in this Letter the quantities defined in 
(§) are at most of the order of 10^^. Also in this range, the error in the derivatives of the 
metric functions is even smaller (the metric functions are very smooth). We also want to 
mention that the Poincare sections shown in this Letter were computed from orbits with an 
accumulated error in the "energy" [cf. Eq. (|^)] smaller than 10^^'^. 

We want to conclude with a discussion of a possible astrophysical implication of our 
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main result: direct rotating particles are less chaotic than counter rotating ones. In a 
rotating center of gravitational attraction we can have structures formed by counter and 
direct rotating particles, thus our result favors larger life times of structures formed with 
direct rotating particles. A deeper discussion of this point will be presented elsewhere. 

The authors thank CNPq and FAPESP for support and to S.R. de Oliveira and N. Santos 
for discussions. 
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FIGURE CAPTIONS 

Fig. 1. Poincare section of test particles moving with angular momentum L = —3.75 in an 
approximate Kerr geometry with rotation parameter a = 0.01 (counter rotation) and mass 
171 = 1. This is a typical section of an integrable system. 

Fig. 2. Poincare section for D = 0.0005, E = 0.965, a = 0, m = 1 and L = ±3.75. The two 
isolated islands around the points (10, 0.05) and (5, 0.075) are parts of the same torus. 

Fig. 3. Poincare section for D = 0.0005, a = 0.01, m = 1 E = 0.965 and L = 3.75 (direct 
rotation). The islands of stability are larger in this case than in the non rotating case (cf. 
Fig. 2). 

Fig. 4. Poincare section with the same parameters of Fig. 3 , except that now L = —3.75 
(counter rotation). The chaotic region increases in a significant way, also the system of 
islands located around the points (10, 0.05) and (5, 0.075) in the two precedent cases has 
disappeared. 
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